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Malaria parasite clones with the highest transmission rates to mosquitoes also tend to induce 
the most severe fitness consequences (or virulence) in mammals. This is in accord with expec- 
tations from the virulence -transmission trade-off hypothesis. However, the mechanisms 
underlying how different clones cause virulence are not well understood. Here, using data 
from eight murine malaria clones, we apply recently developed statistical methods to infer 
differences in clone characteristics, including induction of differing host-mediated changes in 
red blood cell (RBC) supply. Our results indicate that the within-host mechanisms underlying 
similar levels of virulence are variable and that killing of uninfected RBCs by immune effectors 
and/or retention of RBCs in the spleen may ultimately reduce virulence. Furthermore, the 
correlation between clone virulence and the degree of host-induced mortality of uninfected 
RBCs indicates that hosts increasingly restrict their RBC supply with increasing intrinsic viru- 
lence of the clone with which they are infected. Our results demonstrate a role for self-harm in 
self-defence for hosts and highlight the diversity and modes of virulence of malaria. 
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1. INTRODUCTION 

Most pathogens exhibit broad genotypic variation in 
the extent to which they harm the host. Across a 
broad array of pathogen species (bacteria, viruses and 
protozoa), some clones within the species barely affect 
hosts at all, while others cause significant health 
impacts [1-4]. In murine malaria, the variation spans 
the range from mild disease that is quickly cleared by 
the host to severe disease that can kill the host [5]. 
Understanding what shapes this variation in virulence, 
both in terms of the biomedical mechanistic (proxi- 
mate) and the evolutionary (ultimate) determinants, 
is of great public health interest [6]. For example, 
understanding the causes of varied virulence can help 
improve management of the immediate health conse- 
quences of infections caused by different clones and 
may also improve long-term prospects for controlling 
the evolution of virulence [7]. A glossary of terms 
appears in table 1. 

* Author for correspondence (charlotte.metcalf@zoo.ox.ac.uk). 



The standard evolutionary perspective on parasite evol- 
ution is that variation in virulence across clones should be 
shaped by the costs and benefits of virulence for parasite 
transmission [8-10]. The spread and persistence of viru- 
lence mutations may be favoured if they covary with 
transmission of the mutation. However, since host mor- 
tality can interrupt transmission, a range of virulence 
traits may lead to equivalent levels of overall fitness for 
the parasite [11]. While it is broadly agreed that the theor- 
etical framework and empirical evidence on the existence of 
virulence -transmission trade-offs need to be further devel- 
oped and extended [10,12], there are some systems where 
the general predictions appear to be borne out [13-15]. 
In particular, there is considerable evidence linking 
transmission rates to variation in the virulence of murine 
and human malaria clones [9]. In mice, malaria parasites 
that replicate more also produce more transmission 
stages, and cause more severe weight loss and anaemia 
of the host. 

However, despite good evidence that malaria clones 
vary in intrinsic virulence [9] , it is not yet clear whether 
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Table 1. Glossary of terms. 



term 



definition 



bystander killing 



clone 

erythropoeisis 
immunopathology 

immune effectors 



transmission 
virulence 



destruction of uninfected red blood cells 
by the immune system; uninfected 
cells can also be removed from 
circulation by the spleen 

genetically uniform lineage of malaria 
parasites 

production of red blood cells 

destruction of uninfected host tissue by 
the immune system 

cells such as macrophages or proteins 
such as antibodies that are capable of 
damaging both parasites and host 
tissue 

spread of a pathogen from one host to 

another or one cell to another 
degree of harm to the host caused by 

infection; this includes direct damage 

by parasites as well as 

immunopathology 



virulence is achieved in the same way or via different 
mechanisms across clones. For the former, for example, 
all clones might be identical except in terms of the rate 
at which parasites invade uninfected red blood cells 
(RBC), and this rate would determine virulence. For 
the latter, variation might be found in several different 
characteristics of the parasite life cycle that would dif- 
ferentially impact virulence. This distinction has 
interesting implications both in terms of considering 
the impact of treatment and control measures (drug 
treatment, vaccination) within and across clones, but 
also from an ecological and evolutionary perspective. 
The degree of similarity in how parasites are exploiting 
host resources and evading immunity will shape the 
outcome of competitive interactions between clones 
[16,17], as well as the potential for clone coexistence. 
For example, some malaria parasite clones have been 
shown to have a high specificity for particular RBC 
age classes [2], and clones that target the same resources 
(e.g. reticulocytes (young RBCs)) may compete more 
strongly and be more likely to exclude each other than 
clones that show no overlap in resource use; the latter 
may be more likely to facilitate each other [18]. 

There is a further, important layer of complexity. 
Immunopathology has long been recognized in biomedi- 
cine as a key determinant of the degree to which host 
health is affected by parasites, and evolutionary biol- 
ogists are increasingly inclined to agree [10,19-23]. 
Indeed, many of the negative health outcomes that we 
think of as malaria parasite virulence are mediated by 
the hosts' own immune systems [23,24]. For example, 
cerebral malaria represents the combined pathogenesis 
of parasites and inflammation [24,25]. Furthermore, 
the immune system destroys (or removes from circula- 
tion) a substantial number of uninfected RBCs during 
malaria infection [26,27]. However, the destruction of 
uninfected RBCs may actually be contributing to the 
long-term control of the infection [28,29]. Establishing 
whether processes perceived as immunopathology (e.g. 



destruction of uninfected host tissue by the immune 
system) are primarily a component of virulence (e.g. 
collateral damage of parasite replication strategies) 
or primarily a component of defence (e.g. host strate- 
gies for limiting parasite resources) is critical to our 
understanding of host -parasite co-evolution. 

Recently, statistical methods have been developed to 
partition causes of within-host malaria dynamics based 
on regression techniques of current parasitaemia on 
past infected and uninfected RBCs [29] that use 
measures of uninfected and infected RBCs one time- 
step in the past to partition mechanisms determining 
the present. This method enables the study of time- 
varying aspects of both host and parasite contributions 
to parasite dynamics. The methods centre around esti- 
mation of a time-varying quantity analogous to the 
transmission rate in canonical susceptible-infected- 
recovered epidemiological models, known as the effective 
propagation number, P e . Knowledge of P e allows direct 
estimation of the time-varying within-host effective 
reproduction ratio R e [30], or the average number of 
new infected cells per infected cell in a previously 
infected bloodstream. This is a key explanatory variable, 
because R e > 1 indicates that the parasite population is 
growing, whereas R e < 1 indicates that it is shrinking. 
Temporal fluctuations in both the effective propagation 
P e and R e may also reveal changes in the relative impor- 
tance of different control mechanisms through time [29]. 
For example, in an experiment where mice were inocu- 
lated with different parasite densities [31], differences in 
effective propagation around the timing of the peak of 
parasitaemia indicate that adaptive immunity is delayed 
in hosts inoculated with lower parasite doses [29]. Such 
regression-based approaches also allow inference of 
time variation in host regulation of erythropoiesis and 
clearance of uninfected RBCs from circulation by 
immune effectors and the spleen. 

Here, we use these new analytical techniques to explore 
the four following questions: (i) What are the similarities 
and differences among rodent malaria clones in the magni- 
tude and time variation of effective propagation and net 
reproductive number?; (ii) What do these patterns reveal 
in terms of similarity and difference of underlying virulence 
mechanisms across clones?; (iii) Does the extent of immu- 
nopathology vary among clones?; and (iv) Can bystander 
killing, defined as the killing or removal of uninfected 
RBCs by the host during an infection, act as a host defence 
mechanism that ultimately aids control of parasites? 



2. MATERIAL AND METHODS 

2.1. System and experimental set-up 

Laboratory mice frequently serve as model hosts for 
Plasmodium chabaudi chabaudi clones [25]. Every 
24 h, and usually around midnight, infected RBCs 
burst, and the resulting merozoites must locate an 
uninfected RBC within a brief time span. For this 
experiment, eight genetically distinct P. chabaudi cha- 
baudi clones were selected (AD, AJ, AQ, AS, AT, BC, 
CW, ER) based on the range of virulence they induce 
[32]. Infections were initiated with an injection of 10 6 
infected RBCs into five C57BL/6 female mice for each 
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clone. RBC density and proportion infected (from 
which density of infected RBCs can be inferred) were 
measured daily from the fourth to the 20th day post- 
infection. Details are provided in previous descriptions 
of these experiments [33,34]. Figure 1 shows time- 
series of uninfected and infected RBCs for each clone, 
and the relationship between maximum asexual para- 
site density and minimum uninfected RBC numbers. 
Assuming that asexual densities are a good proxy for 
the density of the sexual transmission stages [9], this 
negative relationship indicates that clones that induce 
the highest maximum asexual parasite density also 
induce the most severe anaemia. These relationships 
are consistent and maintained across clones [33]. 
As seen here (figure lc), and in previous work [9], the AS 
and CW clones show substantially less virulence than 
other clones in terms of weight and RBC loss. Materials 
and data are available upon request to C.J.E.M. 

2.2. Parasite effective propagation 

We estimate the effective propagation number, P e , a 
quantity analogous to the transmission coefficient in 
canonical models of between-host infection dynamics 
[29]. Temporal variation in P e is shaped both by the 
effects of RBC age structure (predicted to shape para- 
site burst size, and invasion rates [35,36]) and 
immunity [35,37-40]. By analogy with between-host 
transmission ecology, given the discrete-time nature of 
malaria dynamics within the blood, if I t is the number 
of infected RBCs at time t, and S t is the number of unin- 
fected RBCs at time £, then the expected number of 
infected RBCs observed at time t + 1 will be A t = P et 
S t Jt 5 where P et is the effective propagation number of 
the infection at time t. This number can be thought 
of as the product of contact rates and transmission prob- 
ability, given that a contact has occurred. Taking the log 
on both sides of the earlier- mentioned relationship, we 
can write 

log(E[(7 t+1 )]) = log(P e , t ) + log(J t ) + log(S t ), (2.1) 

where E[(/ t+1 )] denotes the expectation. Applying 
regression techniques to (2.1) enables estimation of P e?t , 
whenever measurements of I t and S t are available. The 
within-host effective reproductive ratio is directly related 
to this quantity via R e t = P e t S t . 

A subset of infected cells may in fact be committed 
to the production of transmission stages, and thus not 
contribute to J t+1 , rather going on to produce gameto- 
cytes. However, we assumed for simplicity that the 
total observed infected cells could be equated with I t 
in (2.1). Although cumulatively, the effect of this allo- 
cation towards sexual stages on the asexual densities 
may be substantial, on the one day time-scale con- 
sidered here, the effect is likely to be very small, 
because the proportion committed is less than a few 
per cent relative to asexual stages [41]. 

2.3. Susceptible red blood cell dynamics 

To also capture RBC replenishment and regulation, we 
can express an equivalent regression model for this 
element of the dynamics. Uninfected RBCs can be 
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Figure 1. Time series of (a) mean uninfected and (b) infected 
RBCs per clone (x 10~ 2 jjlI x ); vertical lines show standard devi- 
ations across five mice; and (c) relationship between maximum 
parasitaemia and depth of the RBC trough (x 10~ 2 jjlI x ) , 
suggesting the virulence -transmission trade-off, with dashed 
lines to indicate quar tiles across individuals. Clone colours on 
the last panel correspond to those on the first two, and range 
from red for the clone resulting in the deepest RBC trough 
(AT), through to deep green for the clone resulting in the 
shallowest RBC trough (CW). 
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lost in three ways: via background mortality, estimated 
at around 0.025 per capita per day [42,43]; via immune 
mediated mortality (i.e. killing of uninfected RBCs 
by the immune system, or retention in the spleen, 
where they can neither carry oxygen for the host, nor 
provide resources for the parasite [26,27]; both of 
these are referred to here as 'bystander killing') and 
by becoming infected. Because we know numbers of 
infected RBCs at every time-point, we can explicitly 
account for this last process, and can write that the 
change in susceptible (uninfected) RBCs, denoted S t , is 

St+i = St + b t — it+i, 

where 7 t+1 is total infected RBCs (encompassing both 
asexual infected RBCs, capable of continued within- 
host infection; and sexual infected RBCs, which 
cannot), and b t is the total change in RBCs attributable 
to processes beyond the parasite (i.e. the balance of ery- 
thropoiesis (generation of new uninfected individuals, 
and mortality, as defined above) . Because we have esti- 
mates for S t and I t at every time- step, b t can be directly 
estimated. Values of b t < 0 indicate that the loss of 
uninfected RBCs due to host-driven processes alone 
exceeds gains from erythropoiesis. Given background 
RBC mortality estimates of around 0.025 per capita 
per day [42,43] in the complete absence of erythropoi- 
esis and taking the starting densities of RBCs (8 x 
10 6 jiT \ fi gure la), the loss of more than 20 uninfected 
RBCs per microlitre per day (8 x 10 4 x 0.025 = 2000 x 
10~ 2 is unlikely without additional mortality 

via bystander killing. We therefore assume that values 
of b t below — 20 [jlI 1 d 1 indicate the occurrence of 
bystander killing. 

To fully resolve how RBC age structure affects sus- 
ceptibility and transmission, studies that distinguish 
reticulocytes from normocytes, and further, to go 
beyond this to estimating the actual age of RBCs in 
days would be invaluable [39]. In the absence of that 
information, we can infer age structure via estimates 
of b t using two boundary assumptions: that parasite 
attack rates are concentrated in either (i) the older, or 
(ii) the younger RBCs. Values of b t > 0 will correspond 
to new RBCs entering into the population. Assuming 
first that infected cell mortality predominantly affects 
older RBCs (normocytes), the values of b t > 0 will 
directly reflect the number of 1-day-old RBCs on day 
t, and the number of 2-day-olds the following day, etc. 
By contrast, for clones that preferentially invade 
younger RBCs, new RBCs entering the population on 
a given day (6 t > 0) are likely to be infected. Thus, to 
reconstruct the age structure, we subtract the number 
of parasites seen on that day from the estimate of b t 
on that day, up to the limit of zero. If there are more 
parasites than the value of 6 t , the remainder can be 
subtracted from the 2-day-olds from the previous 
day, etc. Neither framework will provide an exact reflec- 
tion of the age-structure, because preferences are 
unlikely to be so strict, and for young RBC preferring 
strains, there is no reason to assume 1-day-olds are tar- 
geted first. However, combining both assumptions of 
attack concentrated in either only young or only old 
RBCs, we obtain the full range of potential age 



structures. Given that our overall conclusions are not 
sensitive to either assumption (figure 2 a 7 ), our study is 
robust to this uncertainty. 

2.4- Implications of the host response 
for parasite dynamics 

To explore the implications for every mouse of parti- 
cular profiles of RBC replenishment and destruction 
as defined by 6 t , we can use the parametrized model 
to explore what would have happened had the observa- 
ble uninfected RBC replenishment and depletion 
been modified (e.g. if hosts had been prevented from 
killing uninfected RBCs) by initiating a simula- 
tion with the number of infected and uninfected 
RBCs observed on the earliest day for which both 
measurements are available, and then simulating 
dynamics according to 

^t+i = St - It+i + &*> 

where P et reflects the clone-specific effective propa- 
gation value for time-step t as estimated earlier, and 
6 t * reflects values of b t modified to the net effect of 
bystander killing of uninfected RBCs (by replacing 
negative values of b t by zero). 

The degree to which these simulations capture 
realistic within-mouse malaria dynamics hangs on the 
assumption that fluctuations in P et and b t are indepen- 
dent. If similar processes underlie their dynamics, 
altering b t will result in alterations in P et that the simu- 
lations will not reflect. However, given our limited 
knowledge of how the key processes play out dynami- 
cally, this approach allows us to at least explore the 
'what if scenario of changes in bystander killing, all 
else being equal. 

3. RESULTS 

The overall time-pattern of P e (figure 2a) broadly 
matches that found in a previous analysis focusing on 
the AS clone [29] , with high early effective propagation 
that rapidly declines. However, many clones here also 
show a considerable increase in P e close to the timing 
of the peak of parasitaemia (figure 2a, around day 
10). This is correlated in time with a shift in the age 
structure of RBCs present in the blood towards reticu- 
locytes (figure 2d). The subsequent decline in P e after 
this peak (figure 2a) is likely to reflect engagement of 
adaptive immunity. The effective reproductive ratio 
Re (figure 2c) shows similar patterns, although the 
increase around day 10 is considerably less marked 
because at this point mice are generally anaemic 
(figure la). 

The pattern of uninfected RBC production and 
destruction (figure 26) is also similar to that seen in pre- 
vious analysis [29] with RBC destruction or removal 
dominating for a period during the early exponential 
expansion of the parasite (6 t < 0 across all mice for 
each clone for at least one day before the peak of infec- 
tion) followed by an increase in RBC production (6 t > 0 
from around day nine, presumably in response to 
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Figure 2. For eight clones (legend, colours reflect virulence measured as maximum anaemia), across days post infection (x-axis), 
(a) effective propagation P e (standard errors shown as vertical lines); (b) changes in RBC density (x 10~ 2 jjlI x ) not due to para- 
sites (range across mice in each clone on each day indicated by vertical bars); values less than 0 (horizontal bar) indicate that 
destruction by immune effectors or the spleen exceeds replenishment; and (c) effective reproduction number, R e ; values less 
than 1 (horizontal line) indicate that the parasite population is shrinking; and (d) proportion of reticulocytes at every time 
point obtained from combined RBC and parasite dynamics (see text); either assuming that normocytes are mostly infected 
(solid line) or assuming that reticulocytes are mostly infected (dotted line). 



anaemia) and then generally a return to equilibrium as 
the infection is controlled. Despite these broad simi- 
larities, all clones show variation in timing and 
magnitude of patterns of 6 t , especially after day seven, 
and the response to clones A J and BC deviate quite 
considerably from the others. The differences are signifi- 
cant: likelihood ratio tests for models of 6 t , including a 
random effect to capture mouse to mouse variation, 
and accounting for temporal autocorrelation (fitted 
with Ime in the nlme package in R) reveal a highly sig- 
nificant interaction between 'day' and 'clone' 
(likelihood ratio test: * 2 = 2410, d.f. = 119, p < 

0. 001). For some clones, as many uninfected cells may 
be lost in one time step to bystander killing as are lost 
to parasites at the maximum peak of parasitaemia 
(e.g. for BC and CW, max(— 6 t )/max(/ t ) > 1, figure 3). 

Simulations that do not include bystander killing 
(figure 4a, b) can result in increased virulence in most 
mice, as measured by the depth of the RBC trough 
(figure 4c), but this result in only seen consistently in 
three of the eight clones. This may be partly due to co- 
variation between the magnitude of bystander killing 
and clone maximum effective propagation: the lower 
quant ile of RBC loss across mice is minimal for clones 
with the maximum effective propagation (figure 5a), 

1. e. highest levels of bystander killing are associated 
with mice infected by clones with the highest maximum 
propagation numbers (table 1). 
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Figure 3. The ratio between minimum b t and maximum I t for 
mice within the eight clones. Values greater than 1 (horizontal 
line) indicate that more cells were killed by immunity or reten- 
tion in the spleen in one time-step than by the parasite in one 
time-step, comparing the maximum of each for each mouse. 
The degree to which this occurs varies across clones. 

To facilitate comparison of the effects of bystander 
killing across clones controlling for individual variability 
and clone specific effects on bystander killing, we also ran 
simulations taking median values of b t across all individ- 
uals and clones, using clone-specific time- varying P e 
values, and initiating simulations with the observed 
median starting values of infected and uninfected 
RBCs calculated across all individuals for each clone 
(taking the first day post infection for which both were 
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Figure 4. (Opposite.) The effect of removing bystander killing 
or spleen retention on parasite virulence, (a) Simulated unin- 
fected RBCs and (b) corresponding simulated infected RBC 
numbers (as in figure 1) obtained by taking individual-specific 
starting RBC densities and individual b t values, and setting b t 
to zero for time-steps when b t < 0. (c) Observed minimum 
uninfected RBCs (xl0~ 2 |xl \ x-axis) versus simulated mini- 
mum uninfected RBC obtained as described earlier (?/-axis); 
vertical and horizontal lines indicate quartiles across individ- 
uals; the trough is deeper when no bystander killing is 
implemented for mice infected by several of the clones; excep- 
tions include many mice in the AD and AS clones. 



available, here day four). Comparing median simulations 
with the observed time series indicates that, under 
median values of 6 t , mice infected with virulent clones 
would experience deeper RBC troughs than observed, 
and mice infected with less virulent clones would experi- 
ence shallower troughs than observed (figure 56, blue 
lines and black lines with squares intersect). Since b t is 
broadly inversely correlated with the magnitude of 
clone effective propagation (figures 2b and 5a), the 
median b t will represent an increase relative to the 
observed b t for virulent clones and a reduction for aviru- 
lent clones. This pattern indicates that upregulation of 
erythropoiesis (and/or reduction of bystander killing) 
may counterintuitively worsen mouse health via 
increased anaemia that is caused by parasite-mediated 
RBC lysis. Removing bystander killing from the simu- 
lations resulted in consistently worsened anaemia for 
all clones except AS (figure 56, triangles relative to 
squares, note that AJ also shows only a very slight 
improvement with anaemia). Although there is variabil- 
ity in the magnitude of this effect, it highlights the way 
in which temporal differences in P e alone affect the 
impact of bystander killing on host health. 



4. DISCUSSION 

Our analysis returned two major conclusions. First, 
although bystander killing might initially be considered 
a component of parasite virulence (e.g. with different 
parasite clones triggering different degrees of indiscrimi- 
nate killing of uninfected RBCs by the immune system 
[44]), our results show that in standardized simulations, 
bystander killing actually reduces the maximum degree 
of anaemia experienced by hosts in all but one clone and 
is thus likely to be a component of the hosts' malaria 
control strategy. Second, clone-specific temporal pat- 
terns of effective propagation through time (figure 2a) 
indicate the existence of intrinsic differences in clone 
exploitation characteristics that are not clearly organized 
along the virulence continuum, and that potentially 
reflect some form of niche partitioning [17,35,36]. Inter- 
estingly, however, although day-to-day fluctuations do 
not show a consistent pattern among clones (with, e.g. 
all the most virulent peaking on day 10), the maxima 
do (figure 5a). 

Differential niche partitioning in malaria parasites 
must rely on variation across three main characteristics 
of their life cycle: (i) clones may differ in their rates of 
uninfected RBC invasion and infected RBC burst size 
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Figure 5. Strain differences in patterns of bystander killing and their impact (a) upper quantile of P e for each clone, versus lower 
quantile of bystander killing (corresponding to the greatest loss) across individuals for that clone (xl0~ 2 |xl \ as in figure 26), 
indicating that clones that have highest effective propagation also experience the greatest magnitude of bystander killing (n = 8, 
y = —2846 + 17 8346 255a;, p < 0.05, r 2 = 0.56); (b) observed (blue circles) and simulated minimum uninfected RBCs 
(xl0~ 2 ixl 1 , 2/-axis) for each clone, taking median starting densities across all mice for each clone, and introducing either the 
median b t observed across all mice at that time step (black squares), or the maximum between this median value and 0 
(black triangles, b t > 0). Removing bystander killing increases anaemia in all cases, except for the AS clone, with differential 
magnitude effects across different clones. 



(which may also vary through time); (ii) clones may 
differ in the degree to which the age of RBCs affects 
these two processes; and (iii) clones may differ in the 
degree to which they elicit immune responses, as well 
as the type of immune response they elicit. Unfortu- 
nately, data are not available to allow us to directly 
estimate the effects of RBC age on effective propagation 
[29], because these will be confounded with the impact 
of immunity on parasites (e.g. effective propagation 
might be reduced by immune-mediated elimination of 
infected cells or merozoites, or by changes in the age- 
structure of RBCs towards RBCs with lower burst 
sizes or invasion rates; and these processes will be indis- 
tinguishable). However, we can dissect the temporal 
patterns of P e and b t to identify which components 
are likely to be implicated. 

Even from the very earliest days post infection, when 
differences in immunity among hosts infected with 
different clones might be expected to be negligible, 
and the RBC age structure should be essentially domi- 
nated by older RBCs (because RBC daily mortality 
rates are less than 0.05 per cent [42,43]), the various 
clones show different magnitudes of P e and sub- 
sequently R e . For example, BC, one of the most 
virulent clones (figure la), shows rather low effective 
propagation on the first day of estimation (figure 2a), 
significantly lower than A J, another virulent clone. 
These early time-step differences may partly reflect 
cross-clone variation in invasion rate or burst size in 
normocytes [29] but cross-clone variation in the degree 
to which immune responses are elicited (or, conversely, 
evaded) is likely to be the key; especially, given evidence 
that for doses of a comparable magnitude to those 
introduced here, both innate and adaptive immunity 
are active during the fourth-fifth day [29]. Previous 



work suggests that BC generally induces a high innate 
response relative to other clones (at least A J, AS and 
CW), as measured by interferon (IFN)^ and tumour 
necrosis factor (TNFa) [33], lending further weight to 
the possibility that early differences in P e are due to 
differential immune response induction. 

The large peak in effective propagation occurring 
around day 10 for many clones is a surprising feature. 
It might be due to changes in the immunological con- 
text of the infection with a dip in immune efficacy 
around day 10, mirroring the double peak of immune 
efficacy used in other models [40]. If this is not the 
case, the pattern suggests a role for differential para- 
site responses to differently aged RBCs [36,40], either 
in terms of numbers of merozoites produced after an 
infected RBC bursts, or rates of uninfected RBC inva- 
sion. The proportion of reticulocytes in the blood 
stream around day 10 is likely to be high (figure 2a 7 ) 
because new RBCs will be entering the bloodstream 
as a result of considerably increased erythropoiesis for 
all clones around this time (figure 26, b t > 0 across 
clones for days after day 10). The fact that the peak 
in P e is timed around day 10 suggests either higher 
burst size in reticulocytes, as supported by some exper- 
imental data [40], or higher invasion rates into 
reticulocytes. The opposite was suggested by previous 
models [36], but not borne out by subsequent exper- 
iments, which found them to be about equal [40]. The 
third possibility is that a combination of RBC age 
effects and waning of immune efficacy [40] is responsible 
for the observed peak in P e . 

The importance of age preferences in determining 
malaria clone virulence has been highlighted previously 
[9]. For example, a single mutation that expands the 
age range attacked by Plasmodium berghei yoeli shifts 
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the species from being a mild to a hypervirulent infection 
[45] . However, here our results suggest that links between 
the targeted age range and degree of virulence may be 
rather complex, because many of the avirulent clones 
as well as the virulent clones show an increase in effective 
propagation at around the time that the age structure is 
expected to transfer to being predominantly composed of 
young RBCs (figure 2 a 7 , day 10). By contrast, AS and 
AJ, which show less sensitivity to this transition, and 
therefore might be expected to have lower effective 
propagation (via invasion rate or burst size) in young 
RBCs are at opposite ends of the virulence spectrum. 

Although complex immune patterns that covary with 
RBC age structure cannot be rejected [39,40], patterns 
throughout the time-series are consistent with persistent 
RBC age effects through time, and of relevance to overall 
virulence phenotypes. Clones that have highest P e in 
older RBCs (e.g. AS, AJ, as evidenced by high early 
propagation, and a reduced response to the influx of reti- 
culocytes around day 10, figure 2), reap the benefits of 
extremely high early i? e , with A J in fact producing a 
number of new infected RBCs close to the burst size 
(merozoite number per infected cell) reported for this 
clone [46]. By contrast, clones that apparently favour 
reproduction in younger RBCs do not experience large 
R e values at the time-point corresponding to their largest 
P e , presumably because many other factors are in play, 
such as low RBC densities and strong immunity. 
Again, this does not translate into the virulence profile 
in any straightforward way. 

In the face of this diversity of clone dynamics 
(figure 2a), it is of interest to ask whether the host 
should be employing a single generic defence strategy, 
or a variety of clone-specific strategies. We do not 
have direct information on how host immune systems 
are negotiating this challenge in terms of killing of 
infected cells (e.g. via clone-specific immune killing 
surfaces as in [29]). However, direct inference of ery- 
thropoiesis and bystander killing rates is possible 
(figure 2b). As mentioned earlier, this would indi- 
cate that hosts employ a very similar approach across 
clones in which they reduce the availability of uninfec- 
ted RBCs via bystander killing during the exponential 
growth phase of the parasite (figure 16), and then 
increase erythropoiesis during peak anaemia before 
returning to equilibrium. Beyond this general consist- 
ency, however, there is also variability linked to clone 
identity. Clone characteristics are correlated with the 
magnitude of host response, such that the highest 
levels of bystander killing were identified in mice 
infected by clones with highest rates of parasite 
growth (figure 5a). This specificity of the host's ery- 
thropoiesis and bystander killing response to different 
clones is also borne out by predictable changes in 
trough depth from simulations of 'the median mouse' 
for each clone, with either the same b t imposed on all 
clones for each time-step (set to the time-varying 
median across the entire dataset), or the clone-specific 
median b t (figure 5b). More virulent clones had shal- 
lower troughs in simulations using the across-clone 
median b t value, as this was larger than their clone- 
specific b t (figure 2b)] less virulent clones experience 
the opposite (figure 5b). 



Overall, these patterns suggest that the best strategy 
from the point of view of the host is to downregulate ery- 
thropoiesis (and/or engage in bystander killing) in 
proportion to the magnitude of RBC loss to parasitaemia 
early on. A proximate cue for such a strategy may arise 
from the production of TNFa by cells of the innate 
immune system upon binding malaria waste products 
as well as intracellular debris [47,48]. The data resolution 
available does not allow us to address this in detail, but 
we might expect reticulocyte-preferring parasite clones to 
be particularly vulnerable to reduction of erythropoiesis 
[28], while normocyte-preferring clones should be more 
vulnerable to bystander killing; an interesting avenue 
for future research. 

While our analysis goes beyond early investigations 
that suggested that common mechanisms were regulat- 
ing parasite growth across clones at all phases of the 
infection [9], there are still a number of key questions. 
In particular, our alteration of R e by setting bystander 
killing to zero ignores many subtleties of RBC supply 
and regulation as well as the potential role of age effects. 
Time-points where b t > 0 may also exhibit bystander 
killing or clearance from circulation by the spleen, but 
this signal may be swamped by the degree of erythropoi- 
esis, which will also shift the age structure markedly 
towards younger RBCs, with potential implications for 
clones' effective propagation. We also assume that mor- 
tality of infected RBCs between the time of RBC 
invasion and the census time in the early morning is neg- 
ligible. This is likely to be a robust assumption, because 
early during the cycle, infected RBCs are antigenically 
and structurally similar to uninfected RBCs. It is only 
over the course of their development that changes in 
membrane topology and display of antigens may occur; 
however, further empirical assessment of this and other 
assumptions will be valuable. 

To conclude, for malaria, the virulence -transmission 
relationship is, on the face of it, rather simple — one 
uninfected RBC is lost every time a parasitized RBC 
is gained; so a negative correlation between peak para- 
site density and peak anaemia is to be expected in 
immune competent hosts. However, our analysis indi- 
cates that the underlying mechanisms and the way in 
which different clones line up along this continuum 
are actually rather complex. Effective propagation 
varies in magnitude across clones, and is higher at 
different times for different clones, owing to differential 
RBC age effects and immune killing of infected cells. 
Furthermore, host responses to different malaria 
clones differ in terms of RBC depletion (figure 5a), 
suggesting that clone-specific responses to parasites by 
the host can play a complex role in contributing to 
short-term measures of virulence (i.e. by increasing 
early anaemia) but thereby mediating the overall 
degree of virulence (including cumulative anaemia) by 
reducing R e (figure 4). While dynamical analysis 
enables us to extend early work that focused solely on 
the timing and magnitude of peak of parasitaemia 
[49-51], a key direction for future work is to link the 
broad patterns we have described here to mechanistic 
differences among clones that account for their different 
immunological vulnerabilities and incorporate details on 
how they exploit specific RBC age structures and the 
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host's clone-specific responses. Experiments incorporat- 
ing a broad range of parasite genetic variation in 
inducing immunological responses will be a key part of 
this endeavour [52] , as will reagents developed to manip- 
ulate the immune responses of laboratory mice. Finally, 
our results highlight that a complete understanding of 
natural selection on immunopathology [20] will require 
consideration of the role of immune-mediated killing of 
uninfected tissue in host defence. 
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